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DEFINITION OF PARAMETERS 

Ij = Spacecraft Transverse Moment of Inertia (1 9.57 slug-ft^ ) 

= Spacecraft Polar Moment of Inertia (21.19 slug-ft“ ) 

= Momentum Wheel Polar Moment of Inertia (0.007488 slug-ft^ ) 

= Momentum Wheel Transverse Moment of Inertia (0.003744 slug-ft^ ) 

Byf = Momentum Wheel Damping Constant (0.005 ft-lb-s/rad) 
cOy = Natural Frequency of the Wheel Inertia Rim (98 rad/s) 

Bjj = Damper Damping Constant (7.38 X 10“^ ft-lb-s/rad) 

= Momentum Wheel Spring Constant (52.82 ft-lb/rad) 

= Damper Spring Constant (4.5 X 10"^ ft-lb/rad) 

/ = Displacement Between Damper and Center of Mass Along the z-Axis ( 1 . 1 48 ft) 

m = Damper End Mass (0.01475 slugs) 

= Displacement Between Damper Hinge Point and Center of Mass Along a Transverse 
Axis (0.0820 ft) 

/•j = Length of Damper Pendulum (0.629 ft) 

= Rate of Energy Dissipation in the Momentum Wheel 
= Rate of Energy Dissipation on the Main Spacecraft Body 
cOj = Relative Angular Rate of Momentum Wheel (209 rad/s) 

“1 

= Spacecraft Body Angular Rates 

"zl 

d = Nutation Angle (Angle between the spacecraft polar axis and the spacecraft momen- 
tum vector) 

X = Spacecraft Nutational Frequency 

H = Momentum Wheel Angular Momentum 
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ON THE ATTITUDE STABILITY OF THE SAS A SPACECRAFT 

by 

G. Fleisher and A. K. Sen 
Goddard Space Flight Center 


INTRODUCTION 


The Small Astronomy Satellite (SAS) is a dual-spin system being designed and developed 
for Goddard Space Flight Center by the Johns Hopkins Applied Physics Laboratory. Basi- 
cally, the system consists of a high speed rotor (momentum wheel) housed inside the main 
spacecraft body. The momentum wheel speed is maintained constant by the use of a syn- 
chronous hysteresis motor. The body is spun up about its longitudinal axis and the spin 
rate (1/12 rpm) is maintained by the use of magnetic torque coils. A nutation damper is 
designed intentionally to be an energy dissipator; however, if additional energy dissipation 
exists due to nonrigid elements in the momentum wheel, the attitude stability of the entire 
spacecraft becomes questionable. The purpose of this analysis is, therefore, to identify and 
model energy dissipation for SAS A in the body as well as in the momentum wheel and then 
to analyze the effect of energy dissipation on system stability. 


STABILITY CRITERIA 

In Reference 1 , Likins has shown that the necessary and sufficient conditions for dual- 
spin stability are expressed by 




< 0 , 


( 1 ) 


where and are energy dissipation rates for the main body (due to the on-board 
damper) and for the momentum wheel, respectively. The frequencies X are expressed as 
follows: 
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Because 7^ and are negative quantities, and X,^ must both be positive for abso- 
lute stability under all conditions. It can be seen that if either X^ or X,^ is negative (Xj^ is 
negative for SAS A) and the corresponding energy dissipation rate is large enough, Equation 
1 might exceed zero. If this condition should occur, the nutational motion would build up, 
and, eventually, the spacecraft would tumble. For the SAS A design, X^ is positive; however 
Xj^ is negative. The following paragraphs will show how expressions for and Tj^ are de- 
rived for SAS A. In reality, the rate of energy dissipation for the momentum wheel will 
affect the rate of energy dissipation for the body damper and vice versa. In this analysis, 
and Tjj will be computed independently of each other and will then be combined, with 
the use of Equation 1 , to give the system stability. This is valid, as it can be shown that even 
for an assumed nutation angle of 10 deg, the maximum flexing angle of the inertia ring 
(Figure 1 ) is at least an order of magnitude smaller than the nutation angle. 



Figure 1— Motion of incremental mass on the momentum wheel rim 
due to flexibility of wheel structure. 
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ENERGY DISSIPATION FROM THE MOMENTUM WHEEL (T^^) 


The only source of energy dissipation for the momentum wheel to be considered here 
is the flexibility of the inertia ring support. Because of the structural damping in the ihertia 
ring support, energy will be dissipated as the support both flexes and rotates without ball- 
bearing deformation. If the flexing angle a is assumed to be small, this source of energy 
dissipation can be modeled as a linear two-degree-of-freedom damping system which is 
subjected to an inertial torque. The equations of motion for this energy dissipation from 
the rotating wheel can be written as (see Appendix A) 

+ [K^ + tan 0 sin (cOj - \)t (3) 

+ lK^^, + = -X//tan0cos(Wj - \)t, (4) 

where the constants and K^^f , which refer to the damping (rate) constant and the spring 
constant of the flexible rotating wheel, respectively, were determined experimentally. 

The necessary experiment was performed by the Johns Hopkins University Applied 
Physics Laboratory. The values of B^ and were determined for both air and vacuum 
conditions. The values used in this analysis were determined in air and represent worst-case 
values (maximum destabilizing effect). With the wheel stationary and with the spin axis 
vertical, an impulse was applied to the outer rim, or inertia ring. Instrumentation was set 
up to provide a record of the vertical displacement of the inertia rim as a function of time. 
The motion of the ring indicated the response was essentially that of a second-order system. 
With this assumption, the natural frequency is found from the response curve, and 
is computed: 

(5a) 

Also, from the response curve, is determined, by noting the time at which the initial 
displacement has been reduced by a factor of 1/e; then B^^, is computed as 

B„ = (5b) 

In terms of the constants and , the steady-state solution of Equation 3 can be ex- 
pressed as follows: 


«x = 


XH tan 6 

+ 2Iwt,03,X)^ + Bloi ] ) 


sin (co / - \0). 
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( 6 ) 


From Equation 3, the average energy dissipation rate due to a torque about the x'-axis can 
be computed as 
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(7) 
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tan^e 


2[(K^^ + 2/^^co»2 + 52 oj2] 
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Similarly, the average energy dissipation rate due to a torque about the /-axis can 
be derived from Equation 4. Thus, the total average energy dissipation rate for the momen- 
tum wheel is finally obtained: 


~ ^Wx ^Wy 


= -ATj tan^ 0, 


( 8 ) 


where 


‘ " [(,K„ + 2/„,.w,X)2 + Bioi^] 


ENERGY DISSIPATION CAUSED BY THE PENDULOUS DAMPER (T^ ) 

As in the previous case, with an initial nutation angle of amplitude 0 for the spacecraft, 
the equation of motion for the pendulous damper can be written as (see Appendix A) 


lrw% 


(9) 


The steady-state solution of Equation 9 is given by 


P = - 


tan 0 


IriiKo - /^X2)2 + z?2x2]i/2 


sin (Kt - y). 


' (10) 


Whence, the average energy dissipation rate due to the pendulous damper can be computed 
as 

Tjd = ”^£>4vg 


= -ATj tan^ 0, 


( 11 ) 
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where 


' 2I^[iKu - 4„X2)2 + /i2x2j‘ 

NUMERICAL RESULTS AND CONCLUSIONS 


As a stability criterion, substitution of Equations 8 and 1 1 into Equation I now yields 


—L - — < 0 


( 12 ) 


Equation 2b indicates that the parameter X|^ for the SAS A is negative for all wheel 
speeds of > 0.007 rpm. For this condition, the stability criterion (Equation 12) reduces 
to 





(13) 


Figure 2 sliows the quantiUes K^ and K 2 l^p as functions of the niomentuin wheel 
speed cOj, This plot indicates that, although the SAS A spacecraft will not be unstable for 
any wheel speed in the range 200-2000 rpm,* the degree of stability will be different at 
different speeds. A convenient measure of this degree of spacecraft stability may be ob- 
tained in terms of a safety factor which is computed by taking the ratio of the maximum 
allowable (K^ = 0.0452 for SAS A) to the actual value of at the operating wheel 
speed selected. Thus, at the chosen operating speed of 2000 rpm for the SAS A momentum 
wheel, the value of is found to be 15.'’8 X 10"'* ft-lb/s, and the safety factor is about 30. 

Thus, it can be concluded that there is no potential problem of instability in the SAS A 
due to the flexing of the momentum wheel inertia ring. 


REFERENCES 

1. Likins, P. W. “Attitude Stability for Dual-Spin Spacecraft”, J. Spacecraft ami Rockets 

4(1 2): 1638-1 643, December 1967. 

2. Bainum, P. M., Fuechsel, P. G., and MacKison, D. L. “On the Motion and Stability of 

a Dual-Spin Satellite with Nutation Damping”, AIAA Guidance, Control and Flight- 
Mechanics Conference, Held in Princeton, New Jersey, August 1969, Paper No. 
69-857, American Institute of Aeronautics and Astronautics, New York, N.Y. 

"'In actual practice, the spacecraft may be shown to be stable for the entire range of wheel speeds between 0 and 2000 rpm. 
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Appendix A 


Equations of Motion 

Energy Dissipation by the Momentum Wheel 

Consider an elementary mass dm on the rim of the momentum wheel (Figure 1) which 
is subjected to inertial forces. If i'j', and k' represent unit veetors for the momentum wheel 
with 0 as the center of mass for the whole spacecraft, so that 

r = iVcos0 + jVsin0 + k'z, (Al) 

then the and z' components of the inertial forces may be found, using small-angle approx- 
imations, as 


dF^ = -dm[{(ji^ + cOj)^rsin0] 
rfFj = dm[z + <I)^/'sin0- co^r cos^ 

+ + o)j)rcos0 

+ coj,(co^ -f- wp/-sin0], (A2) 

where dm = phrdrd<f>, p being the mass per unit area of a wheel rim of any arbitrary 
thickness h. 

In order to find the moment about the x'-axis due to the above forces, note that 
z = (r sin 0)o!^ ; thus, the total moment about the x'-axis due to the forces acting on all the 
elementary masses on the wheel rim is given by 

= JJ (r sin (p)dFj - jjzdp 2 

= (2phj^ r^d^ sin^ <l)d(^[\ + 

+ (^Jq sin0cos0f/^ [-Oiy ■+ 


= + (W, + + W(,(C0^ + COj)], 


's' X 


y ' z 


(A3) 
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where = /»(</) + di )/4, m is the mass of tlie wheel rim and r/j and r/j are, respectively, 
the inner and outer radius of the rim. 

Tlie torque given by Equation A3 must equal the sum of the damping and the restoring 
(spring) torques arising from the values of and for the momentum wheel. Thus, 

or, in terms of the angular rates and of the spacecraft body, 

= j [cb^ - (o)^ + 2Wj)Wj^] sin Wj/ + (cb^^ + cosco^/ j 

- 2Wj<Oj^)sinWj/ + (w^ + 2WjC0^) cos co^r), (AS) 

(since co^). 

Now if, as an approximation, the spacecraft is assumed to be freely nutating with a constant 
amplitude 0 (though, in reality, 9 will be a slowly-varying time function because of the 
presence of dissipation), we have 

0 )^ = cos X/ 

<j3y = < 0 ^ sin \t, (A6) 

where 

and H (= is the z component of the total angular momentum vector. 

Thus, by the use of Equation A6, the torque equation given by Equation AS can be 
reduced to 

^WT^x Avr^i^x 

= //tanO(X - 2w )sin(co - X)/ = X// tan 0 sin (w^ - X)/, (A7) 

tj 

since co^ ^ X < co^, 

(/, - 1 ^) 01 , + ^ ^ 

* ■ /y ~ h ' 

and, for the momentum wheel, 

I — 0’ 

‘W ~ -^WT' 
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^Similarly, if the moments about the jy'-axis (for which z = a^r cos0) are considered 
and Equation A2 is used, it can be shown that 

tan 0 cos (cOj - \)t. (A8) 

Pendulous Damper 

For the case of a single pendulous damper, the equation of motion as derived by 
Bainum et al. (Reference 2) can be expressed as 

ij + bJ + = M(w^ - (A9) 

where = w/-j , n = mr^ 1 and P is the damper angle. It is assumed that (inlM) < 1 , s 0, 
and s 0. 

Again, with the angular rates of the spacecraft body as given by Equation A6, the 
equation of motion for the pendulous damper can be reduced to 

U + V + V = 


T 


-ju(X + sin X/ 

u\H tan 0 . ^ . 

- - — 7 sm \ t, 


(AlO) 


(since co^ <X). 
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